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Plasmids have a key role in the horizontal transfer of genes among bacteria. Although plasmids are 
catalysts for bacterial evolution, it is challenging to understand how they can persist in bacterial 
populations over the long term because of the burden they impose on their hosts (the 'plasmid 
paradox'). This paradox is especially perplexing in the case of 'small' plasmids, which are unable to 
self-transfer by conjugation. Here, for the first time, we investigate how interactions between 
co-infecting plasmids influence plasmid persistence. Using an experimental model system based on 
interactions between a diverse assemblage of 'large' plasmids and a single small plasmid, pNI105, in 
the pathogenic bacterium Pseudomonas aeruginosa, we demonstrate that positive epistasis 
minimizes the cost associated with carrying multiple plasmids over the short term and increases 
the stability of the small plasmid over a longer time scale. In support of these experimental data, 
bioinformatic analysis showed that associations between small and large plasmids are more 
common than would be expected owing to chance alone across a range of families of bacteria; more 
generally, we find that co-infection with multiple plasmids is more common than would be expected 
owing to chance across a wide range of bacterial phyla. Collectively, these results suggest that 
positive epistasis promotes plasmid stability in bacterial populations. These findings pave the way 
for future mechanistic studies aimed at elucidating the molecular mechanisms of plasmid-plasmid 
interaction, and evolutionary studies aimed at understanding how the coevolution of plasmids 
drives the spread of plasmid-encoded traits. 
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Introduction 

Plasmids have a key role in bacterial adaptation to 
novel environments by providing bacterial popula- 
tions with access to a pool of accessory genes that 
allow them to adapt to specific environmental 
stressors (Gogarten and Townsend, 2005; Smets 
and Barkay, 2005). For example, plasmids have 
had a pivotal role in the evolution of antibiotic 
resistance in pathogenic bacteria by allowing resis- 
tance genes to rapidly spread between bacterial 
strains (Alekshun and Levy, 2007; Martinez et aL, 
2009). In spite of the key role that plasmids have in 
bacterial evolvability, understanding how they can 
persist in bacterial populations in the absence of 
positive selection is challenging because of the 
following reasons: (i) plasmids generally impose a 
fitness burden on their hosts, producing a fitness 
advantage for plasmid-free segregants (Bouma and 
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Lenski, 1988; Dahlberg and Chao, 2003; Enne et aL, 
2004; Dionisio et aL, 2005; Yates et aL, 2006; San 
Millan et aL, 2010; Humphrey et aL, 2012); 
(ii) sustained selection for plasmid-encoded traits 
can result in the movement of beneficial genes to the 
chromosome, rendering the plasmid redundant 
(Modi et aL, 1992; Bergstrom et aL, 2000; Toleman 
and Walsh, 2011); and (iii) about half of the 
plasmids are incapable of horizontal transmission 
by conjugation (Smillie et aL, 2010). These asym- 
metries are thought to prevent the maintenance of 
plasmids as bacterial parasites in the absence of 
positive selection for plasmid-encoded genes 
(Gordon, 1992; Levin 1993). Hence, the long-term 
persistence of plasmids in bacterial populations is 
known as the 'plasmid paradox' (Harrison and 
Brockhurst, 2012). In this paper, we show that 
positive epistasis between co-infecting plasmids 
minimizes the cost of plasmid carriage and increases 
the ability of plasmids to persist in the absence of 
selection for plasmid-encoded traits, and we argue 
that epistasis may have an important role in 
resolving the 'plasmid paradox'. 

Plasmids can be divided into two groups according 
to their size and biology: 'small' and Targe' 
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(Smillie et aL, 2010). Large plasmids may be capable 
of horizontally transferring themselves via conjuga- 
tion, can potentially carry many beneficial genes 
and use active partitioning systems and post- 
segregational killing systems to prevent segrega- 
tional loss of plasmids. Segregational loss would 
potentially occur at a relatively high frequency 
owing to the low copy number of large plasmids 
(typically ^2). In contrast, small plasmids lack 
conjugative machinery, carry few potentially bene- 
ficial genes and rely on a high copy number 
(typically >10) to avoid segregational loss. There- 
fore, it is especially challenging to understand how 
small plasmids persist in bacterial populations 
because they generate a burden that is similar to 
that produced by large plasmids (San Millan et aL, 
2010), but lack some of the mechanisms for stability 
within a host that large plasmids enjoy. 

Our understanding of the evolution of bacteria- 
plasmid associations is largely based on the analysis 
of selection in systems containing a bacterial host 
and a single plasmid (Stewart and Levin, 1977; 
Freter et aL, 1983; Simonsen, 1991; Bergstrom et aL, 
2000; Lili et aL, 2007; Ponciano et aL, 2007; Svara 
and Rankin, 2011; Harrison and Brockhurst, 2012). 
Co-infection of bacterial strains with multiple 
plasmids is common, but the influence of inter- 
actions between plasmids in evolution remains 
poorly understood (Paulsson, 2002). It is well 
known that the coexistence between plasmids that 
are members of the same incompatibility group is 
unstable because they share common elements 
involved in replication and partitioning (Novick, 
1987). In this case, interactions between co-infecting 
plasmids compromise plasmid stability. However, it 
is also conceivable that interactions between com- 
patible plasmids promote plasmid maintenance by 
either increasing the benefits associated with plas- 
mid carriage or by decreasing its costs. These types 
of positive interactions that ameliorate bacterial 
fitness have been recently shown between chromo- 
somal mutations and conjugative plasmids confer- 
ring antibiotic resistance (Silva et aL, 2011). 

In this work, we investigate the evolutionary 
consequences of interactions that occur between 
co-infecting 'large' and 'small' plasmids using an 
experimental system consisting of the bacterium P. 
aeruginosa and six plasmids that replicate in 
Pseudomonadaceae, including five large plasmids 
and a single small plasmid. We tested how inter- 
actions between large and small plasmids influence 
the costs (decreased competitive ability in the 
absence of antimicrobials) and the benefits (resis- 
tance to antibiotics) of plasmid carriage. To comple- 
ment these short-term competition experiments, we 
carried out a long-term plasmid stability experiment 
in which we allowed populations of plasmid- 
carrying bacteria to evolve in the absence of 
antibiotics for 100 generations to determine how 
interactions between plasmids influence plasmid 
persistence in the absence of selection for plasmid- 



encoded genes. Finally, we analyzed data on the 
distribution of plasmids across all 1866 bacterial 
strains currently available in GenBank to test 
whether patterns of plasmid carriage in nature are 
consistent with the results of our experiments. 

Materials and methods 

A detailed version of the methods including the 
bioinformatic analysis is available as Supplementary 
Material. 



Bacterial strains, plasmids and antibiotic susceptibility 
testing 

The bacterial strains and plasmids used in this 
study are detailed in Table 1. Minimal inhibitory 
concentrations were performed in triplicate accord- 
ing to Clinical and Laboratory Standards Institute 
(2009). 



Competitive fitness assays 

The fitness of each strain carrying the different 
plasmid combinations was determined relative to 
the PAOl-green fluorescent protein (GFP) strain. 
Precultures of the strains were incubated at 37 °C 
with 200 RPM shaking overnight in 3 ml of LB broth 
(Fisher Scientific, Bridgewater, NJ, USA). Precul- 
tures were diluted 100-fold in 3 ml of fresh LB broth 
and incubated in the same conditions until they 
reached mid-exponential phase (OD 600 of ^0.5). 
Cultures of the strains were then mixed at a ratio of 
^50% PAOl/plasmid to 50% PAOl-GFP. The exact 
initial proportions were confirmed via flow cyto- 
metry using an Accuri C6 Flow Cytometer Instru- 
ment (BD Accuri, San Jose, CA, USA). Mixtures 
were diluted 200-fold in fresh LB and competed for 
16 h at 37 °C with no agitation (^7 generations). 
Again, the final proportion was measured by 
flow cytometry. The fitness of the strain carrying 
the plasmid(s) relative to the PAOl-GFP strain 
was determined using the formula 

^plasmid = log e (Nfinal,plasmid/ Mnitial,plasmid) /l°g e (^final,GFp/ Mnitial,GFp) > 

where Wpi asm id is the fitness of the strain carrying the 
plasmid, N initia i ;plasmid and N fina i ; pi asmid are the 
numbers of cells of the strain carrying the plasmid 
before and after the competition and N initia i GFP and 
A/fmai,GFP are the numbers of cells of PAOl-GFP 
before and after the competition. 



Quantification of plasmid copy numbers 
The copy numbers of plasmids were determined by 
quantitative polymerase chain reaction (qPCR) using 
an ABI Prism 7000 sequence detection system (Life 
Technologies, Paisley, UK). In Supplementary Table 
S2, we present the primers and conditions used for 
this analysis. To calculate the copy number of 
plasmid per chromosome, we used the formula: 
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Table 1 Strains and plasmids used in this study 



Strain or plasmid Description and characteristics Reference, Accession no: 



Plasmids 


Group h 


Size 


Resistance 


Transmission 


Origin 


Year 






(kb) 


phenotype c 








pNI105 


NA 


5 


Neo r , Kan r 


Non-transmissible d 


Soil 


1991 


(Itoh et al, 1991) AB032347 


pBS228 


IncP-la 


89 


Str r , Car r , Tet r 


Mobilizable 


Waste water 


1981 


(Haines et al, 2007) 














NC_008357 


Rmsl49 


IncP-6 


57 


Str r , Gen r , Car r 


Mobilizable 


Clinical 


1975 


(Haines et al., 2005) 
















NC_007100 


pAKDl 


IncP-ip 


58 


Str r , Hg r 


Conjugative 


Soil 


1998 


(Sen et al, 2011) JN106164 


PAMBL-1 


ND 


30 


Gen r , Car r , Caz r , 


Conjugative 


Clinical 


2005-2007 


(Tato et al, 2010) 








Mer r , Hg r 








GQ422829 


PAMBL-2 


ND 


30 


Str r , Gen r , Car r , 


Conjugative 


Clinical 


2005-2007 


(Tato et al, 2010) 



Caz r , Mer r GQ422828 
Strains 



P. aeruginosa PAOl 

P. aeruginosa PAOl-GFP 



PAOl/ 
PAOl/ 
PAOl/ 
PAOl/ 
PAOl/ 
PAOl/ 
PAOl/ 
PAOl/ 
PAOl/ 
PAOl/ 
PAOl/ 



pNI105 
pBS228 
Rmsl49 
pAKDl 
PAMBL-1 
PAMBL-2 
pNI105 pBS228 
pNI105 Rmsl49 
pNI105 pAKDl 
pNI105 PAMBL-1 
pNI105 PAMBL-2 



Susceptible strain used as a recipient for plasmids 
PAOl tagged with the green fluorescent protein (GFP) 
gene integrated at the Tn7 insertion site. Gm r , Chl r 
Neo r , Kan r 
Str r , Car r , Tet r 
Str r , Gen r , Car r 
Str r , Hg r 

Gen r , Car r , Caz r , Mer r , Hg r 

Str r , Gen r , Car r , Caz r , Mer r 

Neo r , Kan r , Str r , Car r , Tet r 

Neo r , Kan r , Gen r , Str r , Car r 

Neo r , Kan r , Str r , Hg r 

Neo r , Kan r > Gen r , Car r , Caz r , Mer r , Hg r 

Neo r , Kan r , Str r , Gen r , Car r , Caz r , Mer r 



NC_002516 

(Choi and Schweizer, 2006) 

This work 
This work 
This work 
This work 
This work 
This work 
This work 
This work 
This work 
This work 
This work 



Abbreviations: NA, not applicable; ND, not determined. 

Sequences present in GenBank for PAMBL-1, PAMBL-2 and pNI105 are partial. 

b Small plasmid pNI105 does not have a defined group, it codes for a rolling circle replication-like rep-protein. 

c Gen, gentamycin; Chi, chloramphenicol; Neo, neomycin; Kan, kanamycin; Str, streptomycin; Car, carbenicillin; Tet, tetracycline; Hg, mercury; 
Caz, ceftazidime; Mer, meropenem. r , high-level resistance phenotype. 

d The small plasmid pNI105 showed a conjugation frequency <10" 7 transconjugants per donor colony-forming unit (no transconjugant detected 
in the conjugation assay) when coexisting with the conjugative plasmids pAKDl, PAMBL-1 and PAMBL-2. 



cn = (l+£ c ) ctc /(l+£p) Ctp xS c /S p , where cn is 
the plasmid copy number per chromosome, S c and 
S p are the sizes of the chromosomal and plasmid 
amplicons (in bp), E c and E p are the efficiencies of 
the chromosomal and plasmid qPCRs (relative to 1) 
and Ctc and Ctp are the threshold cycles of the 
chromosomal and plasmid reactions, respectively. 



Results and discussion 

To investigate interactions between plasmids, we 
used six plasmids that are able to replicate in 
Pseudomonas and were recovered from different 
sources, including five large plasmids from a diverse 
range of incompatibility groups: pBS228, Rmsl49, 
pAKDl, PAMBL-1 and PAMBL-2, and a single small 
plasmid pNI105 (Itoh et al, 1991, 2003; Haines 
et al, 2005, 2007; Tato et al, 2010; Sen et al, 2011) 
(Table 1). The five large plasmids range in size from 
30 to 90 kb and confer resistance to a wide variety of 
antibiotics and the heavy metal mercury (Table 1). 
Three of these large plasmids are conjugative 
(pAKDl, PAMBL-1 and PAMBL-2) and two are 



mobilizable (pBS228 and Rmsl49, that is, they are 
not self-transmissible but they can be transferred by 
a trans-acting conjugative element present in the 
same host bacteria). The small plasmid pNI105 is 
only 5kb in length and confers resistance to 
neomycin and kanamycin; this plasmid cannot be 
mobilized by any of the large plasmids used in this 
study under our experimental conditions (see 
Materials and methods). To study interactions 
between large and small plasmids, P. aeruginosa 
PAOl was transformed with the six different 
plasmids independently and with all five small- 
large plasmid combinations (Table 1). 



Positive epistasis minimizes the burden of plasmid 
carriage 

To determine how epistasis between plasmids 
influences the burden associated with plasmid 
carriage, we directly competed plasmid- 
bearing strains against a plasmid-free strain 
(PAOl) in culture medium lacking antibiotics. 
We distinguished between plasmid-bearing and 
plasmid-free cells by using a plasmid-free strain 
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carrying a GFP gene integrated at the Tn7 insertion 
site (Choi and Schweizer, 2006). As a control for the 
costs and benefits associated with this construct, we 
also competed the GFP-tagged strain (PAOl-GFP) 
against an unmarked PAOl strain (the average fitness 
of PAOl-GFP relative to PAOl was 0.977, s.e. = 0.005; 
n = 5). There was significant variation in the cost of 
carriage among large plasmids (one-way analysis of 
variance: P= 0.0002, F = 17.70, df =4, 10), with some 
plasmids lacking a detectable cost and others impos- 
ing a substantial burden of up to 14.3% (Figure la). 
Consistent with previous work, we found that the 
small plasmid imposed a large cost (mean fitness 
cost=13.2%; s.e. = 2.5%; 77 = 3 biological replicates 
of 5 technical replicates each). 

Epistasis between plasmids (e) was calculated as 
previously described by Silva et aL, 2011, using the 
formula e W s .l - W s _ x W_ L , where W s _, W L and 
Ws.l are the fitnesses (relative to the plasmid-free 
PAOl) of the strains carrying a small plasmid, a 
large plasmid and both small and large plasmids, 
respectively, and the error of e [o e ] was estimated by 
the method of error propagation (Trindade et aL, 
2009; Hall and MacLean, 2011; Silva et aL, 2011). 
As in Silva et ah (2011), we interpreted a positive value 
for e — o e as significant positive epistasis between 
plasmids. On the whole, we found evidence that 
strong positive epistasis minimizes the burden 
associated with carrying multiple plasmids. We 
detected significant positive epistasis in 4/5 plasmid 
associations (pNI105 and Rmsl49 e = 0.119 ± 0.065, 
pBS228 e = 0.125±0.05, pAKDl e = 0.159 ± 0.068, 
PAMBL-2 e = 0.085 ±0.076) and in the remaining 
case the interaction between pNI105 and PAMBL-1 
did not differ significantly from the multiplicative 
model [e = 0.005 ±0.068). In the four associations 
showing positive interactions, strains carrying a 
large plasmid paid very low costs or no additional 



costs at all for acquiring the small plasmid, in spite 
of the large cost associated with introducing a small 
plasmid into a plasmid-free strain (Figure lb). 

Co-infection does not reduce plasmid copy number 
Previous work has shown that the burden associated 
with plasmid carriage correlates to plasmid copy 
number (Patnaik, 2000; Harrison et aL, 2012). To test 
the hypothesis that a reduction in plasmid copy 
number drives positive epistasis between plasmids, 
we measured the number of copies of each plasmid 
per cell using qPCR (Figure 2). As expected, the 
copy number of the small plasmid (18 copies/cell; 
s.e. = 2.4) was substantially higher than the copy 
number of the large plasmids, which ranges between 
^1 and 3. Across all of the plasmid-plasmid 
associations, we found that acquiring a large 
plasmid did not significantly change the copy 
number of the small plasmid (two-sample Mests: 
false discovery rate-corrected P>0.4 for all combi- 
nations, see the legend of Figure 2; Benjamini and 
Hochberg, 1995). Similarly, acquiring a small plas- 
mid did not significantly change the average copy 
number of any of the large plasmids (Two-sample t- 
tests: false discovery rate-corrected P>0.1 for all 
combinations, see the legend of Figure 2; Benjamini 
and Hochberg, 1995). In summary, the coexistence of 
multiple plasmids in the same cell did not lead to 
any consistent changes in plasmid copy number, 
demonstrating that reduced plasmid copy number 
cannot explain the low cost associated with carrying 
multiple plasmids. 

Plasmids contribute independently to antibiotic 
resistance 

Recent studies have shown that the cross talk 
between large plasmids and the bacterial chromo- 
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Figure 1 Fitness of plasmid-carrying strains, (a) Fitness (±s.e.) of PAOl carrying one and two plasmids, relative to the plasmid-free 
ancestor, P. aeruginosa PAOl. The horizontal dotted red line represents the fitness of PAOl. Vertical dotted black lines separate each pair 
of strains carrying a large plasmid alone and combined with the small plasmid pNI105. (b) Positive epistasis between large and small 
plasmids. Comparison of the relative fitness of P. aeruginosa PAOl carrying large and both large and small plasmids. The blue line 
represents the expected fitness if the fitness cost produced by large and small plasmids is multiplicative. The red line represents the 
expected fitness if the small plasmid (pNI105) produces no extra cost. Error bars represent standard error. The points above the blue line 
represent associations of plasmids showing positive epistasis. 
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I 
I 



• pNllOS 

• pBS228 

• Rmsl49 
pAKOl 

o PAMBL-1 
PAMBL-2 



PA01/ PA01/ PA01/ PA01/ PA01/ PA01/ PA01/ PA01/ PA01/ PA01/ PA01/ 
pNI105 pBS228 pNI105 Rms149 pNI105 pAKD1 pNI105 PAMBL-1 pNI105 PAMBL-2 pNI105 
pBS228 Rms149 pAKD1 PAMBL-1 PAMBL-2 

Figure 2 Plasmid copy numbers in PAOl-transformed strains. Average number of plasmid copies per cell ( ± s.e.). Copy numbers were 
determined by qPCR of plasmid-encoded genes using the chromosomal gene rpoD as a reporter. The copy number of the small plasmid is 
not altered by the coexistence with a large plasmid (two-sample f-tests: pNI105 and Rmsl49: P= 0.785; pNI105 and pBS228: P= 0.454; 
pNI105 and pAKDl: P= 0.454; PAMBL-1: P = 0.935; PAMBL-2: P= 0.684; df = 10). The presence of the small plasmid does not change the 
copy number of large plasmids (two-sample f-tests: Rmsl49: P= 0.104; pBS228: P= 0.104; pAKDl: P = 0.271; PAMBL-1: P= 0.932; 
PAMBL-2: P = 0.272; df=4). 



some produces alterations in gene expression that 
minimize the fitness cost associated with plasmid 
carriage (Harr and Schlotterer, 2006; Navarre et ah, 
2006; Doyle et al, 2007; Shintani et al, 2010; 
Yun et ah, 2010). An alternative explanation for 
the reduction in the cost of coexisting plasmids that 
we observed could be that plasmid gene expression 
decreases during co-infection. To test this hypo- 
thesis, we measured the resistance of strains carrying 
single plasmids or pairs of plasmids to antibiotics 
and the heavy metal mercury, as defined by their 
minimal inhibitory concentration (Supplementary 
Figure Si; Supplementary Table Si). All of the 
plasmids in this study carry antibiotic resistance 
genes (Table 1) that degrade or inactivate antibiotics, 
and it is well established that phenotypic levels of 
resistance for this class of resistance genes are well 
correlated with enzyme expression levels (Uhlin 
and Nordstrom, 1977; Mick et al, 2008). We found 
that coexisting plasmids contributed additively to 
antibiotic resistance, as demonstrated by the fact 
that acquiring a plasmid that did not confer 
resistance to a particular antibiotic did not change 
the level of resistance (that is, minimal inhibitory 
concentration) of cells harboring a plasmid confer- 
ring resistance to that antibiotic (paired £-test: 
P=0.752, t= -0.318, df=44, Supplementary 
Figure Si). When both plasmids conferred resis- 
tance (pNI105 confers high-level resistance to 
kanamycin and neomycin, and PAMBL-1 and 
PAMBL-2 confer very low-level resistance to kana- 
mycin and to kanamycin and neomycin, respec- 
tively), the strains carrying both of them showed an 
additive effect on their resistance level 
(Supplementary Table Si, Supplementary Figure 
Si). These results provide good evidence that 
positive epistasis between plasmids does not arise 
as a consequence of the repression of plasmid- 
encoded antibiotic resistance genes. In the broader 



context of antibiotic resistance, this implies that 
acquiring multiple plasmids provides bacteria with 
a mechanism to become resistant to multiple 
antibiotics without incurring a large fitness cost. 
Coexistence of antibiotic resistance plasmids produ- 
cing positive epistasis might therefore have a key 
role in the maintenance of multidrug resistance 
phenotypes in nature. These results could help 
explain why multidrug resistance in clinical patho- 
gens is sometimes associated with carrying multiple 
resistance plasmids (San Millan et al., 2009; Granier 
et al, 2011; Wei et al, 2011). 



Positive epistasis favors plasmid stability 
As plasmids impose substantial fitness costs, it is 
challenging to explain how plasmids can be stably 
maintained over long time scales. In theory, positive 
epistasis between plasmids could help maintain the 
stability of plasmids in the absence of selection for 
plasmid-encoded genes by minimizing the benefits 
associated with losing individual plasmids. Speci- 
fically, under positive epistasis, co-infecting plas- 
mids mask each other's costs so that the loss of 
either plasmid individually will generate only a 
marginal fitness benefit to the bacterium. Conver- 
sely, if plasmids do not interact (multiplicative 
effects on fitness), plasmid stability would not be 
increased because the benefit associated with losing 
a plasmid is the same irrespective of whether or not 
the plasmid co-occurs with a second plasmid 
(Figure 3). 

To test the hypothesis that plasmid coexistence 
increases stability, we analyzed the stability of the 
plasmids both alone and in combination. We 
allowed strains of bacteria carrying pNI105 either 
individually or in combination with one large 
plasmid to propagate in an antibiotic-free culture 
medium. After 100 generations of growth in the 
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Figure 3 Positive epistasis increases plasmid stability. A sche- 
matic representation of the fitness of a bacterial strain carrying 
zero, one or two plasmids. In this example, the strain carrying two 
plasmids can present either multiplicative fitness (no epistasis 
between plasmids, e = 0) or positive epistasis between plasmids 
(e>0). The arrows represent the change in fitness associated with 
the loss of a plasmid. The strain carrying two plasmids that 
produce multiplicative costs (e = 0) will experience the same 
fitness advantage associated with the loss of each of the plasmids 
as in the strain that carries the single plasmids (orange and blue 
arrows, the increase in fitness is the same). On the other hand, in 
the case of positive epistasis, the advantage associated with the 
loss of a plasmid will be lower than in the previous cases (orange 
dotted arrow); the loss of the plasmid would even produce a 
decrease in fitness if the fitness of the strain with two plasmids is 
higher than the fitness of the strain with only one of the plasmids 
(blue dotted arrow). 
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Figure 4 Stability of the small plasmid, pNI105, after evolution. 
Stability of pNI105 (percentage of plasmid-bearing cells, ± s.e.) in 
the populations with the different plasmid combinations after 100 
generations of evolution. PAOl/pNI105 is the control strain, 
carrying pNI105 alone. The plasmids of the PAOl/pNI105 
PAMBL-1 strain display multiplicative fitness effects. The 
remaining plasmid combinations display positive epistasis and 
are grouped together. In the four populations where pNI105 is 
both associated with a large plasmid and also displays positive 
epistasis, there are no significant differences in the stability of 
pM05 among them (one-way analysis of variance, P= 0.193, 
F = 1.73, df=3, 20). 'Average' is the average of the results (±s.e.) 
of the stability of pNI105 in the strains showing positive epistasis 
between plasmids. 



absence of antibiotics, we measured the percentage of 
plasmid-bearing cells in each population (for every 
plasmid) as a measure of plasmid stability (Figure 4). 
Large plasmid stability was high. Rmsl49, pAKDl, 
PAMBL-1 and PAMBL-2 were present in ^99% of the 
colonies tested both when they were alone and in 
combination with pNI105. pBS228 showed a lower 
stability (80.3% alone and 87% in combination with 
pNI105) probably owing to the particularly low copy 
number per cell of this plasmid in PAOl (Figure 2). 
The stability of the small plasmid was high, but 
variable, and, as expected, we found differences in 
the stability of pNI105 across treatments (analysis of 
variance: P= 0.0185, F = 3.645, df=2, 33). The 
stability of pNI105 was lowest (mean = 91.7%; 
s.e. = 1.2%; 77 = 6) in the control treatment that only 
carried pNI105. In the four strains carrying combina- 
tions of plasmids that displayed positive epistasis, the 
average stability of pNI105 (95.9%; s.e. = 0.93%; 
27 = 4) was significantly higher than in the pNI105- 
alone control strain (Dunnett's test: P= 0.016, 
t= 2.493, df=2, 33), consistent with the hypothesis 
that positive epistasis between plasmids increases 
plasmid stability. Further evidence to support this 
hypothesis comes from the observation that the 
stability of pNI105 when coexisting with PAMBL-1, 
a plasmid that showed almost perfect multiplicative 
fitness with pNI105 (mean stability = 93.3%; 
s.e. = 1.98%), is not significantly different from the 
stability in the pNI105-alone control strain 
(Dunnett's test: P= 0.325, t= 0.773, df=2, 33). 

In general, the results of the stability experiments 
showed an overall tendency for higher stability for 
plasmid pNI105 when coexisting with large plasmids 
with which it exhibits positive epistasis. Taking into 
account that pNI105 is not mobilizable in this model 
system, these results support the hypothesis that 
positive epistasis may have a role in increasing 
plasmid persistence in bacterial populations. 



Small and large plasmids are associated in nature 
One of the limitations of this work is that we only had 
access to one small plasmid conferring resistance to 
antibiotics for the PAOl model. Therefore, we decided 
to extrapolate beyond our experimental system using a 
bioinformatic approach to test for signatures of 
increased plasmid stability stemming from positive 
epistasis between large and small plasmids. We 
analyzed the patterns of distribution of plasmids across 
the genomes of free-living bacteria available in 
GenBank. If positive epistasis between large and small 
plasmids has a key role in the maintenance of small 
plasmids, then we expect the following: (a) strains 
carrying only small plasmids should be under-repre- 
sented, (b) strains carrying a combination of large and 
small plasmids should be over-represented and/or they 
should carry a disproportionately high frequency of 
small plasmids. 

To analyze the distribution of large and small 
plasmids, it was necessary to establish a size limit 
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between those two groups. In the complete genome 
database (1866 genomes, see methods), this limit is 
not evident (Supplementary Figure S2) owing to the 
heterogeneity in plasmid biology and sizes among 
the different bacterial taxa. Specifically, there is a 
positive correlation between chromosome and plas- 
mid lengths (Smillie et ah, 2010) that masks the 
specific distribution of plasmid size within each 
group when all the groups are pooled. Therefore, we 
analyzed the distribution of size within bacterial 
families where plasmid biology and genome size are 
more homogeneous. Those families having ^20 
genomes carrying plasmids were further analyzed: 
Enterobacteriaceae (174 genomes), Bacillaceae (76 
genomes), Staphylococcaceae (42 genomes) and 
Lactobacillaceae (43 genomes). In three of these 
families (Enterobacteriaceae, Bacillaceae and Sta- 
phylococcaceae), the distribution of plasmid by size 
is bimodal and a clear limit could be established to 
perform the analysis (Figure 5; for the distribution in 
the Lactobacillaceae family, see Supplementary 
Figure S3). The plasmid-bearing genomes of these 
three families were analyzed. Three categories of 
strains were designated: strains with small plasmids 
(one or more, no large plasmids), strains with large 
plasmids (one or more, no small plasmids) and 



strains with both large and small plasmids (one or 
more of each). The expected frequencies of strains in 
each category, as well as the small and large plasmid 
frequencies within each category, were calculated 
for each family under the null hypothesis that 
associations between large and small plasmids are 
random. We then tested our two hypotheses and 
found a deficit in strains carrying only small 
plasmids both in Enterobacteriaceae [% 2 test, 
P< 0.0001, x 2 = 15.80, df=l) and Staphylococca- 
ceae (x 2 test, P= 0.034, x 2 = 4.51, df=l), but not in 
Bacillaceae (x 2 test, P= 0.523, % 2 = 0A1, df=l) 
(Figure 5b). The number of strains carrying both 
large and small plasmids was not significantly 
increased in any of the families (x 2 test: Enterobac- 
teriaceae, P= 0.501, % 2 = 0A5, df=l; Bacillaceae, 
P= 0.709, x 2 = 0-14> df=l; Staphylococcaceae, 
P=0.95, x 2 = 0.01, df=l). On the other hand, the 
frequency of small plasmids in the strains carrying 
both large and small plasmids was higher than 
expected by chance for both Enterobacteriaceae (x 2 
test, P< 0.0001, x 2 = 41.67, df=l) and Bacillaceae 
(X 2 test, P=0.02, x 2 = 5.35, df=l), but not for 
Staphylococcaceae (x 2 test, P=0.25, % 2 = 1.34:, 
df=l) (Figure 5c). Therefore, at least one of 
the two predictions supporting the tendency for 
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small-large plasmid associations was confirmed in 
each of the three different families tested. Collectively, 
these bioinformatic results are consistent with the 
hypothesis that positive epistasis between large and 
small plasmids is common, and that interactions with 
large plasmids could have a key role in stabilizing 
small plasmids. This analysis reveals that in three of 
the most clinically relevant bacterial families there is 
a tendency for associations between small and 
large plasmids, especially in Enterobacteriaceae. 
Interestingly, these families present numerous 
examples of small plasmids carrying antibiotic 
resistance genes (Bernhard et al, 1978; Cao et al, 
2002; Bjorland et al, 2003; Hauschild et al, 2006; 
Zioga et al, 2009; Chen et al, 2010; Kadlec et al, 
2010; Wang et al, 2012), and thus these associations 
could lead to the acquisition of multiresistance 
phenotypes, which represent a serious concern for 
infectious disease therapeutics. 

Alternatively, it could be argued that this tendency 
for associations between large and small plasmids 
arises from the fact that mobilizable small plasmids 
cannot transmit horizontally unless they co-occur 
with a large plasmid. Previous analysis based on the 
presence of mobilization (relaxases) genes showed 
that approximately half of the small plasmids are 
supposed to be non-mobilizable (Smillie et al, 2010), 
and thus they would not be able to co-transfer with 
large plasmids. Using the same approach, we found 
that exactly 50% (70) of the small plasmids in these 
three bacterial families are non-mobilizable, and the 
distribution of non-mobilizable small plasmids does 
not differ significantly (x 2 test, P= 0.783, x 2 = 0.08, 
df = 1) between genomes with only small plasmids (8 
out of 17) and genomes with both large and small 
plasmids (62 out of 123). This analysis suggests that 
the tendency for the coexistence of small and large 
plasmids is, at minimum, not solely driven by 
increased horizontal transmission of small plasmids 
in the presence of large plasmids. The caveat to this 
analysis is that defining the potential for plasmid 
mobilization based on relaxases genes may under- 
estimate the mobilizable capacity of small plasmids 
as they could be using trans- acting or as-yet 
unknown relaxases. 



Co-infection with multiple plasmids is common 
Although our experimental work focused on study- 
ing interactions between one small and five large 
plasmids, it is conceivable that positive interactions 
shape the evolution of compatible plasmid coex- 
istence regardless of size. If positive epistasis 
between plasmids is a general phenomenon, then 
we expect that the frequency of strains carrying a 
single plasmid will be under-represented relative to 
strains carrying multiple plasmids. 

To test this prediction, we analyzed the number of 
plasmids present in the 1866 bacterial genomes that 
are currently available in GenBank. Our null 
hypothesis is that if plasmids are distributed 



randomly across bacteria, then the number of 
plasmids per genome will be Poisson distributed. 
Our justification for this null hypothesis is that if n 
observations (in this case, plasmids) are randomly 
distributed across x discrete intervals (in this case, 
genomes), then the distribution of observations per 
interval (in this case, plasmids/genome) will follow 
the Poisson distribution (Zar, 1974). The Poisson 
distribution provides a useful statistical null model, 
because a significant deviation from the Poisson 
distribution implies that chance alone cannot 
explain the observed distribution of plasmids across 
bacterial genomes and that other forces, such as 
selection and conjugation, must influence this 
distribution. We found that the distribution of 
plasmids per strain deviated significantly from the 
Poisson expectation (Figure 6a, x 2 test* P< 0.0001, 
X 2 = 1045.31, df = 4), which is perhaps not surpris- 
ing because selection will influence the distribution 
of plasmids. For example, we found that strains that 
lack plasmids are more common than we would 
expect owing to chance (x 2 test, P< 0.0001, 
X 2 = 271.43, df= 1), and strains carrying one plasmid 
are less common than expected (x 2 test, P< 0.0001, 
X 2 = 305.31, df=l), which is consistent with selec- 
tion acting against strains bearing plasmids owing to 
the cost of plasmid carriage. Given that selection 
against plasmids is effective, strains carrying multi- 
ple plasmids should become progressively more 
under-represented as plasmid number increases if 
fitness interactions between plasmids are multi- 
plicative or negative. In contrast to this idea, we 
found that strains carrying multiple (>2) plasmids 
are over-represented (x 2 test, P< 0.0001, x 2 = 136.66, 
df=l); specifically, we found that the overrepre- 
sentation of strains carrying multiple plasmids 
increases exponentially with plasmid number, as 
we would expect if positive epistasis occurs 
between plasmids (analysis of variance: P= 0.0003, 
F = 78.784, df=l, 5, R 2 = 0.928) (Figure 6b). 
Although these findings are consistent with positive 
epistasis facilitating multiple plasmid coexistence, 
it is important to keep in mind that the majority of 
plasmid-positive strains still carry only one or two 
plasmids (454 out of 629). 

It could be argued that the tendency for plasmids 
to coexist in these genomes is a consequence of the 
bias toward sequencing strains of clinical origin, 
where antibiotics impose a strong pressure and 
preferentially select bacteria carrying coexisting 
antibiotic resistance plasmids. To address this 
possibility, we performed the same plasmid dis- 
tribution analysis for all bacterial phyla with ^30 
genomes sequenced and ^10 genomes carrying 
plasmids: Proteobacteria (819 genomes), Bacteroi- 
detes (72 genomes), Actinobacteria (229 genomes), 
Cyanobacteria (48 genomes) and Firmicutes (475 
genomes). These phyla differ in the proportion of 
clinically relevant species that they contain; in fact, 
the phylum Cyanobacteria includes no species of 
clinical significance at all. Regardless of the 
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representation of clinically relevant species, the 
same main features were observed in every group 
when compared with what was predicted by the 
Poisson distribution (Figure 6c): (i) a nonrandom 
distribution of plasmids [% 2 test, P< 0.0005), (ii) a 
higher frequency of strains with zero plasmids [% 2 
test, P<0.001, df=l), (iii) a lower frequency of 
strains with one plasmid [% 2 test, P< 0.0001, df=l) 
and (iv) a higher frequency of genomes with >2 
plasmids [% 2 test, P<0.05, df=l). In Actinobacteria 
and Bacteroidetes, this last analysis could not be 
performed because the expected frequency of strains 
with >2 plasmids was lower than five, but the 
observed number was considerably larger than 
expected in both cases (Actinobacteria: expected = 1.56, 
observed = 7. Bacteroidetes: expected = 1.45, 



observed = 9). Finally, in four out of the five phyla 
(Proteobacteria, Bacteroidetes, Cyanobacteria and 
Firmicutes), strains carrying two plasmids were 
more rare than expected [% 2 test, P<0.02, df=l), 
whereas in Actinobacteria they were more frequent 
than expected [% 2 test, P<0.001, df=l). Notably, in 
all the different phyla, the same general tendency 
toward plasmid coexistence was observed (regard- 
less of plasmid size), indicating that this is a general 
phenomenon occurring across the bacteria phylo- 
geny and independent of the antibiotic pressure 
experienced by bacterial pathogens. Interestingly, 
the distribution of plasmids in bacterial hosts that 
we observed is similar to the classic negative 
binomial distribution observed for metazoan para- 
sites infecting individual hosts in a population 
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(Shaw and Dobson, 1995; Poulin, 2013). Other 
biases may be present in the databases, such as the 
loss of plasmids during genome assembly, and thus 
even when we consider that the bioinformatic 
analyses support the hypothesis that positive 
epistasis has a role in plasmid stability in bacterial 
populations these results should be interpreted with 
caution. 



Conclusion 

It is important to understand the underlying evolu- 
tionary biology of bacteria-plasmid associations 
because of the key role that plasmids have in 
bacterial evolution. The central unresolved problem 
in our comprehension of bacteria-plasmid interac- 
tions is to understand how plasmids can persist in 
bacterial populations in the absence of selection for 
plasmid-encoded genes. Previous models and experi- 
ments have focused on tackling this problem in 
systems containing a bacterium and a single plasmid 
(Stewart and Levin, 1977; Freter et al, 1983; 
Simonsen 1991; Levin 1993; Verger et al, 1993; 
Bergstrom et al, 2000; Lili et al, 2007; Ponciano 
et al, 2007); ironically, we find that strains carrying a 
single plasmid are relatively rare (Figure 6), and that 
the coexistence of multiple plasmids in the same 
strain is relatively common (Figure 6). Our experi- 
mental work in P. aeruginosa shows that positive 
epistasis between plasmids increases plasmid stabi- 
lity, and bioinformatic analysis makes the tantalizing 
suggestion that positive epistasis may have a key 
important role in resolving the 'plasmid paradox' 
across bacteria. We hope that this study will pave the 
way for future mechanistic studies aimed at under- 
standing why positive epistasis occurs between 
plasmids, and evolutionary studies aimed at under- 
standing how positive epistasis impacts the coevolu- 
tion of bacteria and plasmids. 
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